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A novel scheme for the creation of a convergent, or focussing, fast-electron beam generated from 
ultra-high-intensity laser-solid interactions is described. Self-consistent particle-in-cell simulations 
are used to demonstrate the efficacy of this scheme in two dimensions. It is shown that a beam 
of fast-electrons of energy 500 keV - 3 MeV propagates within a solid-density plasma, focussing at 
depth. The depth of focus of the fast-electron beam is controlled via the target dimensions and 
focussing optics. 


The study of fast-electron generation and subsequent 
transport in high-density plasmas is important for nu¬ 
merous applications including proton and ion beam pro¬ 
duction [I], isochoric heating of high density matter for 
opacity and equation-of-state studies [2], x-ray sources, 
and fast-ignition inertial fusion |3]. 

There are a number of schemes to ignite a small mass of 
inertially-confined Deuterium-Tritium (DT) fuel, initiat¬ 
ing a thermonuclear burn wave which propagates through 
a larger mass of fuel. The most developed inertial con¬ 
finement fusion (ICF) scheme is central hotspot ignition, 
in which the spherical implosion of a plastic capsule con¬ 
taining a solid DT layer compresses and heats a small 
mass of DT gas to ^ 5 keV and a density-radius prod¬ 
uct (pr) of 0.4 gem -2 . In order to create these hotspot 
conditions simultaneously via spherical compression, ex¬ 
cess energy must be expended compressing the cold DT 
ice and ablator, limiting the fusion energy gain for more 
robust central hotspot target designs. The principal al¬ 
ternative laser-driven ICF schemes are shock and fast ig¬ 
nition, these may require less energy to form the hotspot 
and hence have the potential for higher gain. 

In the fast-ignition scheme, the DT fuel is firstly com¬ 
pressed using long-pulse lasers, then a shorter duration, 
high-intensity 1 x 10 20 Wcm -2 ) laser is used to gen¬ 
erate a beam of fast-electrons, these stop collisionally 
within the compressed DT ice, heating it. The coupling 
efficiency from the high-intensity laser to the DT fuel is 
determined by the fraction of laser energy absorbed into 
fast-electrons, their energy spectrum, divergence, and the 
distance from the laser absorption surface to the com¬ 
pressed DT target or core [4|. A number of studies [5] 
have found the inherent electron beam divergence to be 
so large that the high-intensity laser energy required for 
fast-ignition is excessive from both a practical and eco¬ 
nomic perspective. This in turn has led to proposals nu 
which seek to address this issue by re-directing the di¬ 
vergent fast-electron beam towards the compressed core, 
thereby reducing the high-intensity laser energy required 
for ignition. 

This letter describes a technique to control the diver¬ 


gence of the fast-electron source. Particle-in-cell simu¬ 
lations are used to demonstrate that this technique can 
be used to create a convergent fast-electron source. The 
fast-electron beam focal length can be controlled simply 
by altering the focussing optics and interaction geometry, 
and it is shown that the generated fast-electron beam can 
be focussed within a solid. 

The Lorentz force equation F = — q(E + v x B) de¬ 
scribes the interaction of a single electron with a plane 
electromagnetic wave, here q is the electron charge, E the 
electric field vector, B the magnetic field vector and v is 
the electron’s velocity. The electromagnetic wave firstly 
accelerates the electron in the direction of the electric 
field, if the magnitude of the electric field is sufficient for 
the electron to obtain a velocity approaching the speed of 
light c, the v x B component of the Lorentz force equation 
becomes significant with respect to that due to E, causing 
the electron to also be accelerated in a direction parallel 
to the wave’s Poynting vector. Defining the normalised 
vector potential as ao = qPi/m e cu: where m e is the elec¬ 
tron mass and u the wave’s angular frequency, we find 
that for a wave with its Poynting vector in the z direc¬ 
tion, and an electric field aligned with the y axis varying 
as E = Eocos(ujt), the y and z components of the elec¬ 
tron’s momentum are respectively: P y = aom e csm(ujt) 
and P z = {a^m e c/A) cos 2(ut). Hence for ao > 2 or 
/A 2 > 5.5 x 10 18 Wpm 2 cm 2 the z component of mo¬ 
mentum exceeds that in y. For even higher intensities, 
the electron’s final trajectory approaches that of the elec¬ 
tromagnetic wave’s Poynting vector, as P z oc ag/4 while 
P y oc ao- This electron acceleration mechanism is well 
known in high-intensity laser plasma physics, where the 
necessary intensities are readily achievable, and is com¬ 
monly referred to as j x B acceleration. 

It follows from the above that for an electromagnetic 
wave which is non-planar, but smoothly varying in the 
y direction (using the above coordinates), the electron’s 
velocity vector will initially be locally parallel to the E 
field. Therefore wherever the electron is when its veloc¬ 
ity reaches ^ c, and hence the v x B component of the 
Lorentz force becomes large, the electron will be accel- 
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Figure 1. The setup used in the PIC modelling and con¬ 
cept illustration. Green depicts the Poynting flux of the fo¬ 
cussing incident laser which enters from the left, red is Au 
density. The solid blue upward arrow indicates the trajec¬ 
tory an electron might initially follow on one cycle of the laser 
when being principally accelerated by the laser’s electric field. 
Wherever the electron is along its upward trajectory when it’s 
speed reaches ~ c, the v X B component of the Lorentz force 
will be perpendicular to both the electric field and the target 
front surface (blue dashed line). In this example the electron 
reaches c at the upward pointing solid blue arrowhead. As the 
radius of curvature of the laser and target front are 20 pm, the 
electron is accelerated toward the focal point at y=0, x=20. 
Conversely on the next half-cycle of the laser the electric field 
accelerates the electron downwards (cyan solid arrow), but 
again the v x B component (cyan dashed arrow) accelerates 
the electron towards the focal point. It should be noted that 
were an electron to reach c at any point along its path over 
the target front surface, the v x B component would always 
be normal to its trajectory at that point. 


erated in a direction tan~ 1 {P z /P y ) = £an -1 (ao/4) with 
respect to the wave’s local Poynting vector, as shown in 
figure [l] For an electromagnetic wave with a phase front 
of uniform curvature (a focussing, spherical wave), the 
Poynting vector at each location on the phase front points 
towards one point - the focal point of the wave/centre of 
curvature. Therefore if ao is sufficient, any electrons ex¬ 
periencing the focussing wave will be accelerated towards 
the wave’s centre of curvature. This geometry occurs 
naturally with a focussing laser beam, prior to the beam 
reaching focus. Hence by combining a focussing laser 
beam with a solid target which has the same radius of 
curvature as that of the laser at the target front, but that 
is convex (as viewed by the laser), it should be possible, 
based on single electron dynamics, to create a conver¬ 
gent or focussing fast-electron source. This interaction 
geometry is illustrated in figure [lj with the Poynting flux 
shown in green while the target density is shown in red. 
Note the target front surface is placed well before the 
laser reaches peak focus, and the laser’s electric field is 
parallel to the y-axis at y=0. 

In order to investigate this interaction geometry, two 
dimensional simulations have been performed with the 
particle-in-cell code EPOCH m • They were initialised 


with 800 particles/species/cell, with convergence tests 
up to 1600. Using cubic spline particle shapes, numer¬ 
ical heating and grid convergence tests (minimum grid 
spacing 1.6 x 10 _3 pm) indicate that a spatial resolu¬ 
tion of up to 150 times the Debye length (0.01pm) is 
acceptable; over 150fs the plasma electrons numerically 
heated from an initial temperature of 89eV to 89.5eV 
(HOeV) with collisions off (on) while the fast-electron 
generation and transport was quantitatively unaffected 
with respect to higher grid resolution runs. The timestep 
was 0.01 fs, while the simulated density was solid density 
Gold/Diamond-like-Carbon/DT with an assumed ionisa¬ 
tion state of 10+/5+/1+. A density scale length of 0.1pm 
was used at the target front surface, varying this altered 
the fast-electron temperature and absorption fraction. 
The laser wavelength was 1.064pm and the peak time- 
averaged Poynting flux was 1 x 10 2 ° Wcm -2 at the tar¬ 
get surface, although the scheme was robust over a range 
of relativistic intensities. Thermal boundary conditions 
were used. In order to prevent anomalous boundary ef¬ 
fects associated with the fast-electron beam, the simula¬ 
tion was stopped before the beam reached the box edge. 

In order to characterise the fast electron source with 
minimal additional complexity, initial EPOCH simula¬ 
tions were performed with all particle collisions off and 
no ion species. As EPOCH updates Maxwell’s equations 
using the current, having no ions is equivalent to hav¬ 
ing infinitely massive ions. Figures [ 5 ] (a)-(d) depict the 
flux of fast-electrons (defined throughout as those with 
kinetic energy greater that 0.5 MeV and less than 3MeV) 
induced by the laser plasma interaction depicted in figure 
[l] (solid Au target). The curved peaks in the flux intensity 
distribution have a spacing of 0.5 pm, or half the laser 
wavelength, indicating absorption by the j x B mecha¬ 
nism. Figures [2] (a)-(d) illustrate the propagation of a 
convergent fast electron beam within the solid density 
target. The radius of curvature of the target front and 
the laser are 20pm, while the fast electron beam reaches 
maximum focus at approximately this depth, as expected 
from the above single-electron dynamics arguments. The 
inset vector plots also confirm the directionality of the 
fast-electrons. Note that in the inset of fig. [2] (a) the ve¬ 
locity vectors outside the solid surface (those to the left 
of the blue dashed line) are not locally perpendicular to 
the target surface. This is because at the time of this 
snapshot they have not yet received the full j x B push 
of this particular cycle of the laser, and hence are still 
predominately being accelerated down the target surface 
by the electric field. 

The physics in the above PIC simulation is simplified 
by the effective use of infinitely massive ions and by omit¬ 
ting collisions. Nevertheless it confirms that in this limit, 
and for the conditions examined, the interaction is rea¬ 
sonably described by single electron dynamics which pre¬ 
dicts both the temporal behaviour and trajectories of the 
fast electrons. Separate simulations (not shown) consis- 
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Figure 2. (a) Fast-electron flux at 18 fs, note the spatial scale is smaller on this plot for clarity, the inset depicts the fast-electron 
flux vectors of a small region (grey box), (b) Fast-electron flux at 50 fs. (c) Fast-electron flux at 75 fs. (d) Fast-electron flux 
at 100 fs, inset depicts the fast-electron flux vectors of a small region (grey box). In all cases the blue dashed line denotes the 
initial position of the solid Au surface. 


tently confirmed that the focal length of the fast-electron 
beam approximately equals the radius of curvature of the 
laser/target front. 

By neglecting ion motion, effects caused by deforma¬ 
tion of the critical surface or expansion due to heating 
of the front surface ions are ignored. Deformation of the 
critical surface may de-optimise the interaction geometry, 
potentially degrading the fast-electron source direction¬ 
ality, reducing the focussed fast-electron beam intensity. 
Furthermore compression or expansion of the front sur¬ 
face is expected to change the fast-electron temperature. 

Figure [3] (a) depicts a simulation with is identical to 
that shown in fig. [2] except that ion motion has now been 
enabled. It can be seen that for the timescales exam¬ 
ined, the fast-electron source is similar to that generated 
when ion motion is prevented. Figure [3] (b) depicts the 
target front surface at 300 fs, from an otherwise identi¬ 
cal run but with reduced target depth for computational 
efficiency. The dense region of the target front surface 
is largely unperturbed, despite the effects of hole-boring 
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Figure 3. (a) Fast-electron beam flux generated from simu- 
lated laser-solid interaction including ion motion but no col¬ 
lisions at 100 fs. (b) Zoomed in density plot of the target 
surface after 300fs. 


and the fact that this interface is Rayleigh-Taylor muni 
unstable due to the opposing gradients in density and 
pressure at the target front surface. This indicates that 
little degradation in the source directionality due to de¬ 
formation of the critical surface caused by ion motion 
would be expected up to this time. The principal differ¬ 
ence observed when ion motion is included is a temporal 
variation in the fast-electron energy distribution function 
(spectrum). This is caused by the laser compressing the 
pre-imposed density ramp, steepening the density profile 
and reducing the absorption into fast-electrons. As the 
front surface ions heat up and expand, a shelf forms in 
the ion density profile m- This causes an associated 
temporal variation in both the fraction of laser-light ab¬ 
sorbed and the fast-electron energy distribution function. 
Despite the additional complexity caused by the ion mo¬ 
tion, this result provides confirmation that a convergent 
fast-electron source is created using this novel interaction 
geometry. 

In order to asses whether a convergent fast-electron 
beam propagates within the target from the convergent 
fast-electron source, the effects of collisions must be in¬ 
cluded during the transport of the fast-electrons through 
the target. Figure [4] (a) depicts the 0.5 — 3 MeV fast- 
electron flux from the same simulation as that shown 
in figs. [5] (e) and [ 3 ] (a) but particle collisions are now 
included. This has been simulated using a Coulomb log¬ 
arithm of both 2 and 4, with little quantative difference 
observed. The fast-electron flux of the former is depicted 
in figure |4^a), again the inset shows the fast-electron flux 
vectors near the source point in approximately the cor¬ 
rect direction for a focussing beam. In comparison to 
the previous runs it can be seen that the extent to which 
the beam converges is reduced, this is perhaps unsurpris¬ 
ing as collisions will cause the beam to ‘bloom’. Figure 
4 (b) shows a lineout along the horizontal axis of fig. 
4^a). The increase in fast-electron flux as a function of 
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Figure 4. (a) Fast-electron beam flux for PIC runs including 
both ion motion and collisions. Inset depicts fast-electron flux 
vectors within the region in the grey box shown in (a), (b) 
Red line is a lineout of fast-electron flux along the y=0 axis. 
Blue line is the same as the red, but smoothed. 


distance into the target indicates the beam is conver¬ 
gent, alt hough this is convolved with temporal variations 
in the laser absorption and the fast-electron energy dis¬ 
tribution function. In this collisional case a greater de¬ 
gree of velocity dispersion of the individual fast-electron 
bunches occurs, so that by 20 pm the 2uo bunching of the 
individual fronts has largely disappeared. As the incident 
laser intensity is 1 x 10 20 Wcm -2 , and the fast-electron 
flux ~ 1 — 3 x 10 19 Wcm 2 the fast-electron absorption 
fraction is of the order of 10-30%. When the Au was tran¬ 
sitioned to a fully ionised DT plasma at target depths 
greater than 1pm, very similar results were obtained. If 
however the front surface of the target was made from 
either Diamond-like-Carbon or DT, deformation due to 
hole-boring caused the fast-electron beam to rapidly de¬ 
grade. 

Based on Bell et al’s resistive self-collimation criterion 
na, it would be expected that the fast-electron beam 
generated using this scheme would be self-confined, how¬ 
ever this is not observed in these simulations. Very 
large magnetic fields 10 4 T) are generated by the 
fast-electron beam, however these comprise many small- 
scale filamentary structures aligned with the local direc¬ 
tion of the fast-electron beam propagation, the size of 
these are not inconsistent with the Weibel instabilityjTSj. 
The small scale magnetic field structures are at least 
in part due statistical noise; even with 1600 parti¬ 
cles/species/cell full convergence of the magnetic field 
was not possible. The RMS amplitude of the magnetic 
field filaments is inversely proportional to particle num¬ 
ber and is described by B Zrms (T ) = —171/n(pc/s) + 1611, 
where pels is the total number of particles per cell. It is 
notable that hybrid-PIC simulations cause these filamen¬ 
tary field structures to largely disappear m, further¬ 
more when smoothed in post processing, the filaments 
merged into a large-scale azimuthal field with the cor¬ 
rect geometry to self-collimate. It should be noted that 
despite the above, quantitative convergence of the fast- 


electron flux was demonstrated with ^ 400 electrons and 
50ions per cell. 

A convergent fast-electron beam has the potential to 
improve the energetics of the fast-ignitior scheme. By re¬ 
ducing the fast-electron beam radius, the hotspot mass is 
reduced and hence the required ignitor beam energy re¬ 
duced. Furthermore by focussing the fast-electron beam 
it may be possible to use a lower laser intensity in order 
to heat the hotspot within the inertial confinement time, 
this would reduce the fast-electron temperature, improv¬ 
ing collisional coupling. Whether the scheme outlined in 
this letter will translate into a viable ignitor beam is as 
yet unknown. This is due to the greatly increased tempo¬ 
ral and spatial scales required to impart sufficient energy 
into the hotspot and simulate the anticipated offset from 
the cone-tip to the compressed DT fuel. In principal this 
fast-electron focussing scheme should be equally applica¬ 
ble to three spatial dimensions, here optimisation might 
be achieved by employing a radially polarised laser beam 
focussing towards a convex target. These issues will be 
addressed in future work. 

In conclusion, a novel scheme for the generation of a 
convergent fast-electron source has been outlined. Proof 
of principal two dimensional PIC modelling has shown 
that this scheme works for the spatial and temporal scales 
examined. Furthermore it is shown that the fast-electron 
beam can propagate and focus within a solid density 
plasma. 
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